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A dark force can impact the cosmological history of dark matter (DM), both explaining observed cores in 
dwarf galaxies and setting the DM relic density through annihilation to dark force bosons. For GeV - TeV DM 
mass, DM self-scattering in dwarf galaxy halos exhibits quantum mechanical resonances, analogous to a Som- 
merfeld enhancement for annihilation. We show that a simple model of DM with a dark force can accommodate 
all astrophysical bounds on self-interactions in halos and explain the observed relic density, through a single 
coupling constant. 



/. Introduction: The paradigm of cold, collisionless dark 
matter (DM) has been extraordinarily successful in explain- 
ing astrophysical observations of structure, from the recom- 
bination epoch to the present large scale structure of the Uni- 
verse. Although all evidence for DM is from its gravitational 
influence, it is expected that DM possesses some type of in- 
teractions beyond gravity. Non-gravitational interactions are 
required to produce DM particles in the early Universe, and 
ultimately determine the DM density observed today. 

Despite its great success, it is unclear whether cold, col- 
lisionless DM can successfully account for the small scale 
structure of the Universe, which may indicate other interac- 
tions besides gravity play a role in structure formation. Pre- 
cision observations of dwarf galaxies by THINGS show DM 
mass distributions with flat cores, compared to steep cusps 
predicted by collisionless DM simulations The gravita- 
tional effect of massive baryonic outflows from supernovae 
can potentially flatten central DM cusps |2), but it is unknown 
whether this effect can explain the observed cores in other 
less luminous (more DM-dominated) dwarf galaxies |3|. An- 
other discrepancy is the apparent underabundance of Milky 
Way (MW) satellite dwarf galaxies, compared to predictions 
from collisionless DM simulations |4 |. The missing low mass 
satellites may simply be fainter than expected if energy injec- 
tion from astrophysical processes strips away interstellar gas 
and suppresses star formation |5|. However, this mechanism 
cannot explain the apparent absence of the most massive sub- 
halos predicted by simulations |6| which are "too big to fail" 
in star formation and are too dense to host any observed MW 
satellite, according to their predicted stellar circular veloci- 
ties |7|. 

These small scale structure anomalies can be explained if 
DM, denoted X, is self-interacting |8|. Recent N-body simu- 
lations have shown that a DM self-interaction cross section 
per unit mass ax/mx ^ 0.1 — lOcm'^/g can flatten the 
central density within dwarf galaxies and subhalos to solve 
the core/cusp problem |9, 10 1. Moreover, the most massive 
subhalos can be reconciled with the observed MW satellites 
since stellar circular velocities are reduced in their central 
cores 121 (To). At the same time, a variety of constraints 
from larger scales (e.g., halo shapes of elliptical galaxy and 
cluster halos ifTTTl . the Bullet cluster 1121 . subhalo evapora- 
tion IIT3I ) have suggested that a-r/mx must be smaller on 



these scales, motivating a velocity-dependent force lfT4l [TSll 
that gives ax/rnx ^ lOcm^/g on dwarf scales |j9|, but is 
suppressed on larger scales. However, simulations with a con- 
stant cross section have shown that the aforementioned con- 
straints are in fact much weaker than previously thought, and 
a constant ar/rrix ^ 0.1 cm^/g is sufficient to solve small 
scale structure anomalies while evading other bounds 1 10|. 

Given these results, it is important to explore the particle 
physics nature of DM self-interactions. For typical weakly- 
interacting DM models, self-scattering has a weak-scale cross 
section ax ^ 10^^^ cm^, far too small to play a role in 
galactic dynamics. An MeV-scale dark force mediator (de- 
noted (p) is needed to give a much larger scattering cross sec- 
tion, (It 1 cm^ {rtix /g) « 2 x 10^^'' cm^ (nix /GeV), re- 
quired to leave observable signatures on DM halos lfT4HT8l . 
A perturbative calculation for ct from <j) exchange gives 
ax ~ 47rQ;^m^/TO^ at small velocity, where ax is the "dark 
fine structure constant," so that 
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in the desired range. However, this calculation breaks down 
for 7710 < axmx, and nonperturbative effects become im- 
portant. These effects have not been studied in general and yet 
are crucial for connecting dark forces to small scale structure. 
In particular, DM self-scattering exhibits quantum mechanical 
resonances, analogous to resonant Sommeifeld enhancement 
for annihilation, as we show below. 

In this Letter, we present a simple model where a dark force 
can simultaneously set the DM abundance and solve small 
scale structure anomalies. In the early Universe, XX — >■ 
provides an efficient annihilation channel for obtaining the 
relic density during freeze-out. This same coupling resolves 
structure problems through scattering on small scales while 
remaining consistent with bounds on MW and cluster scales. 
We consider both symmetric and asymmetric DM models 
which involve attractive and repulsive DM self-interactions. 
In calculating the scattering cross section, we take a numeri- 
cal approach and cover the full parameter space including the 
nonperturbative quantum mechanical regime which has not 
been explored before. We show that resonances can arise for 
a wide range of DM mass, mx ~ GeV — TeV. Furthermore, 
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our numerical calculation confirms analytical formulae widely 
used in literature for computing ctt in the classical and Born 
limits. 

//. DM Annihilation and Elastic Scattering: We consider a 
Dirac fermion DM particle X, coupled to a dark force vector 
boson (p with mass via 

^int = gxXrx^^, (2) 

where gx is the coupling constant. We assume that X is 
weakly coupled to the SM (e.g., through kinetic mixing of 
(j) with U{1)y hypercharge) so that X thermalizes with the 
visible sector in the early Universe f20l. 

DM freeze-out is governed by the velocity-weighted anni- 
hilation cross section for XX (pel), given by (crw) an ~ 
Tra^/ml- where ax = g\/{ATi). For symmetric DM, 
where DM consists of equal densities of X and X, we re- 
quire ((Tu)an ~ 6 X 10^^^ CTo^ / s to obtain the observed relic 
density. For asymmetric DM, the present DM density is deter- 
mined by a primordial asymmetry between X and X, in anal- 
ogy to the baryon asymmetry. In this case, we require larger 
(fw)an to deplete the symmetric X, X density, leaving behind 
only the residual asymmetric X density as DM. Thus, we have 
ctx ?i 4 X 10^^ (mjf/GeV), with the lower bound saturated 
for symmetric DM. Asymmetric DM allows for a broader re- 
gion of parameter space, since annihilation XX — > (fxf) suffi- 
cient to set the reUc density only places a lower bound on a.x, 
rather than fixing it to a particular value as a function of mx ■ 

In our model, the same dark force carrier also medi- 
ates DM self-interactions. Here, the relevant quantity is the 
scattering cross section weighted by the momentum transfer, 
(Jt — J dft (1 — cos 0) da/dil, where da/dO, is the usual dif- 
ferential cross section. The nonrelativistic interaction between 
two DM particles mediated by (p is described by a Yukawa po- 
tential 

F(r) =± — e-™*'-. (3) 
r 

Since (/> is a vector, XX XX scattering is repulsive (+), 
while XX XX is attractive (— ). For symmetric DM, both 
attractive (X-X) and repulsive {X-X or X-X) interactions 
are present; for asymmetric DM, where DM consists of only 
X after the freeze-out, self-interactions are only repulsive. 

Since both scattering and annihilation occur through a com- 
mon interaction, the cross sections are related. When 
is massless, the scattering cross section scales roughly as 
arp ~ {o'v)an/iJ*- If this relation holds to dwarf scales 
(v ^ 10 km/s), the transfer cross section is (jT/mx ^ 
10"^ cm^/g (TeV/mx), which is too large compared to that 
preferred by the simulation results 1 9 10 1 unless the DM mass 
is larger than 100 TeV. Therefore, a nonzero is essential, 
softening the velocity-dependence of at small v due to the 
finite range of the dark force. 

The calculation of for a Yukawa potential with ^ 
is non-trivial. We collect analytical results, where applicable, 
in the appendix. Within the Born approximation (valid for 
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FIG. 1: Velocity-dependence of or for sample parameters within 
different regimes. Blue line shows Bom formula (|4|, in agreement 
with numerical results (blue dots), for mx ~ 4 GeV, — 7.2 
MeV, ax ~ 1-8 x 10^*. Green line shows classical formula js}, in 
agreement with numerical results (stars), for mx = 2 TeV, m^ — 1 
MeV, ax ~ 0.05. Red lines show or in the resonant regime for 
mx ~ 100 GeV, ax ~ 3.4 x 10""^, illustrating s-wave resonance 
(solid, m^ = 205 MeV), p-wave resonance (dot-dashed, m^ = 20 
MeV), and s-wave antiresonance (dashed, m^ — 77 MeV). 



ctxmx /"rriij) ^ 1), ctt can be computed perturbatively. Out- 
side the Born regime, multiple scatterings lead to a nonper- 
turbative modification of the DM two-body wavefunction, and 
an analytical approximation has been obtained only within the 
classical limit {mxv/mtp ^ 1). However, outside the Born 
and classical regimes, no analytic description is possible, and 
one must compute ctt by solving the the Schrodinger equa- 
tion numerically using a partial wave analysis |16, 19 1. Within 
this "resonant" regime, ctt exhibits a rich structure of quan- 
tum mechanical resonances (for the attractive potential case).' 
Computing within this regime is crucial for understanding 
for what parameters a dark force can explain simultaneously 
small scale structure problems and the DM reUc density. 

To illustrate the different regimes and behaviors of DM self- 
scattering. Fig. [T] shows aT/mx as a function of v for an 
attractive potential, for several parameter choices. The blue 
(green) line shows the analytic result for for a parameter 
point within the Born (classical) regime (see appendix); these 
formulae are in excellent agreement with our numerical re- 
sults, shown by the blue dots (green stars). The red lines cor- 
respond to three parameter points within the resonant regime. 
The solid red line shows an s-wave resonance, with ax grow- 
ing as v~'^ at small velocity. The dot-dashed line shows a p- 
wave resonance, where ctt shows a resonant peak at finite v. 



' Ref. 1 16 1 previously studied tliis effect for limited parameter choices moti- 
vated by cosmic ray excesses; here, we have adopted a more efficient nu- 
merical procedure (described in a forthcoming publication 1191 ^ allowing 
us to explore the full parameter range in detail. 
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FIG. 2: Symmetric (left) and asymmetric (rigiit) DM parameter space in mx-m^ plane. Blue regions show where DM self-scattering solves 
dwarf-scale structure anomalies, while red (green) lines show bounds from Milky Way (cluster) scales. Numerical values indicate {ot) /mx 
in cm'^/g on dwarf ("dw"), Milky Way ("MW"), and cluster ("cl") scales. For symmetric DM, ax is fixed to obtain the observed relic density; 
for asymmetric DM, ax = 10~^ is fixed to deplete X, X density for nix ^ 300 GeV (dotted line). Dashed lines show extrapolation using 
analytic formulae, while "x" marks parameter points utilized in Fig. 1. 



These two cases illustrate how cft may be enhanced at dwarf 
scales due to resonances. The dashed line shows an exam- 
ple with an antiresonance (the Ramsauer-Townsend effect), 
which can suppress at small v. All of these parameters 
have been chosen to give the correct DM relic density and 
cr/'mx ^ 0.1 — lOcm^/g to solve structure problems on 
dwarf scales (except for the antiresonance case). 

///. Results: We show the complete parameter space where 
a dark force can account for DM small scale structure and 
relic density. For scattering, to compare with astrophysi- 
cal bounds, we consider the velocity-averaged cross section 
[ot) ~ J d^v (Jt e-^^^ /*'o/(27rw§)'^/^, where wo is the most 
probable velocity for a DM particle. Fig. [2] shows contour 
plots of (ctt) for two cases, symmetric and asymmetric DM, 
in the mx-m^ parameter space. 

For symmetric DM (Fig.|2j /e/f), we take the average of at- 
tractive and repulsive cross sections, <tt = (crf^' + (/r^^)/2, 
with ax chosen to reproduce the observed DM relic density 
at each point. ^ The blue contour regions show {ot) /mx on 
dwarf scales (wq — 10 km/s) in the ranges 0.1 — Icm^/g 
(light) and 1 — lOcm^/g (dark) to solve small scale structure 
problems. The lower range is prefered for a constant cross 



section; Ref. ifTOl found 0.1 cm^/g matched small scale struc- 
ture observations, while 1 ctd? /g caused too low central den- 
sities in dwarf spheroidals. Simulations with a t;-dependent 
classical (attractive-only) force prefered the upper range (or 
larger) [9]. The red (green) contours show (ot) l^x = 0.1 
and 1 cm^/g on MW (cluster) scales with vq = 200 (1000) 
km/s, showing the approximate upper limits from observa- 
tions. Ref. [10| found that Icm^/g produced a too-small 
central DM density in galaxy clusters and is only marginally 
consistent with MW-scale halo shape ellipticity constraints, 
while 0.1 cm^/g is consistent with these constraints |10|. In 
the resonant regime, we have computed ax numerically. This 
region shows a pattern of resonances for mx ^ 10 GeV 
- TeV, where erf." is enhanced, allowing for larger mx for 
fixed {(jT)/mx- The dashed lines indicate where we use 
analytic formulae to extrapolate our results into the Born 
{mx ^ iTL^/ax) and classical (mx ^ m^fi/v) regimes. Our 
numerical calculation maps smoothly into these regions, again 
confirming our agreement with the analytic formulae.^ The 
crosses show the example parameters from Fig.[T]for the res- 
onant (mx = 100 GeV), Born (mx = 4 GeV), and classical 
(mx — 2 TeV) regimes. 

Most of these resonant features correspond to s-wave res- 



^ We compute the relic density by solving numerically the Boltzmann equa- 
tions for DM freeze-out, accounting for a possible Sommerfeld enhance- 
ment in {(jv). We assume X kinetically decouples at a temperature 
0.5 MeV, e.g., if X were weakly coupled to electrons [20). 



' The small discrepancy on cluster scales is because (o"t) at these parame- 
ters is dominated by phase space with v <^ vq, where the classical approx- 
imation is not valid, even though mxvo/ra^ 2> 1- 
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onances, and their location in parameter space is given ana- 
lytically by nix ~ TT^n'^mci,/ (Qax), where n = 1,2,3, etc. 
This condition was derived for Sommerfeld enhancements in 
annihilation II2TI . but the same bound state formation arises in 
scattering as well. Taking ax ~ 4 x 10^^(mx/GeV) to fix 
the relic density, we obtain nix ~ 6 A GeV {m,f, /MeVY ^^n. 
This condition matches the locations of resonances in our nu- 
merical results. 

For asymmetric DM (Fig. |2] right), we take a repulsive- 
only cross section, ctt — cr^^, and no resonances occur. We 
fix ax ~ 10^^, which provides sufficient depletion of the 
symmetric X, X density for mx ^ 300 GeV (dotted line); 
above this line, ADM freeze-out would require an additional 
annihilation channel or a larger ax (which changes the (ctt) 
contours). Numerical and analytic results for (ctt) /mx are 
indicated as in the symmetric case. 

From sub-GeV to multi-TeV DM mass, our results show 
that a dark force can successfully explain both DM struc- 
ture and the DM relic density, for ~ 100 keV — GeV. 
The mx ^ GeV region corresponds to the Born limit; here, 
contours at different vq converge, indicating that ax is ap- 
proximately constant in v. At larger mx, ctt is more sup- 
pressed at larger vq. Therefore, possible evidence for DM 
self-interactions on cluster scales |.22,| may point toward light 
DM. 

IV. Conclusions: Dark forces may play an important role 
in the dynamics of DM, analogous to electromagnetic or nu- 
clear forces in the visible sector. We have shown that a sim- 
ple, generic model with a dark force can simultaneously ex- 
plain the DM relic abundance during freeze-out and solve 
small scale structure anomalies in dwarf galaxies and subha- 
los, while satisfying constraints on larger galaxy and cluster 
scales. We have presented a comprehensive picture of the 
parameter space of our model, considering both symmetric 
and asymmetric DM, with attractive or repulsive dark forces. 
Within the full parameter space spanning these different cases, 
we have shown that the DM relic density and self-scattering 
can accommodate a wide range of DM and mediator masses. 
Importantly for narrowing this range, future astrophysical data 
favoring or more strongly excluding self-interactions on larger 
scales would prefer mx ^ GeV or mx ^ GeV, respectively. 
However, N-body simulations over a larger parameter region, 
including within the resonant regime, would be necessary for 
detailed comparison with observations. Experimental tests 
may also detect the dark force directly, depending on its cou- 
pling to visible matter ||23]| . 
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Appendix: We collect analytic formulae for (Tt used in 
the literature. In the Born regime (axmx /m.^ ^ 1), for 
both attractive and repulsive forces, a perturbative calculation 



gives ifMl 

m^v^ 

where R = mxv/mtp. In the classical regime {mxv/m^ ^ 
1), a solution to the classical equations of motion gives for an 
attractive potential |14, 24| 

r i|/32in(i + r0 p<io-' 

a^J-^J (1 + l-5/3^-'') 10-1 </3< 103 

(5) 

where P = 2axm^/{mxv'^), and for the repulsive case ll25l 

r ^/?2ln(l + r') /3<1 

clas ^ I "^4> /gN 

^ ~ \ ^( In 2/? -In In 2/3)' P>1 ' 
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